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Abstract 

A theory of linearly elastic orthotropic shells is presented, with potential ap- 
plication to the continuous modeling of Carbon NanoTubes. Two relevant fea- 
tures are: the selected type of orthotropic response, which should be suitable 
to capture differences in chirality; the possibility of accounting for thickness 
changes due to changes in inter-wall separation to be expected in multi-wall 
CNTs. A simpler version of the theory is also proposed, in which orthotropy is 
preserved but thickness changes are excluded, intended for possible application 
to single- wall CNTs. Another feature of both versions of the present theory is 
boundary-value problems of torsion, axial traction, uniform inner pressure, and 
rim flexure, can be solved explicitly in closed form. Various directions of ongoing 
further research are indicated. 

Keywords: shell theory, single- and multi-wall carbon nanotubcs, torsion, trac- 
tion, pressure, and rim-flexure problems 



1 Introduction 

The application that motivated this work is the modeling of carbon nanotubes 
(CNTs). When CNTs are employed as nanodevice components, they are re- 
garded as elastic beam-like or shell-like objects and their mechanical response is 
characterized in terms of an as-small-as-possible number of stiffness and inertia 
parameters. To define and evaluate these parameters is the common goal of 
all modelers; a way to achieve it is to try and bridge between the microscopic 
scale of molecular mechanics and the macroscopic scale of continuous structure 
mechanics, by way of a mesoscopic scale, at which concepts from discrete struc- 
ture mechanics apply. At the onset of putting together a bottom-up model of 
this sort [T], we realized that ordinary shell theories, which presume an isotropic 
three-dimensional response of the material comprising the shell, could not pos- 
sibly guarantee an accurate macroscopic account of the mesoscopic texture of 
single- wall carbon nanotubes (SWCNTs): a glance to armchair and zigzag CNTs 
(FigureQ]) suggests instead an orthotropic response in planes orthogonal to radial 
directions (see Figure^ where the three little cylinders suggest what probes one 
should cut out of a cylindrical shell-like body in order to determine its material 
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moduli). Moreover, when modeling multi-wall carbon nanotubes (MWCNTs), 
it seems to us important to allow for thickness distension: we conjecture that 
thickness changes are essentially due to changes in the inter-wall distances, as a 
consequence of the interplay of the applied loads with the van der Waals inter- 
actions between adjacent walls. A search of the literature convinced us that we 
better produced such a shell theory ourselves. This paper describes the results 
of our efforts, results that have been partly anticipated in [4] and that - so we 
believe - may find application also in contexts different from the mechanics of 
nanotubes. 

In the next section, we expound the general lines of a theory of linearly 
elastic orthotropic shells of constant referential thickness 2e, whose geometry 
(Section I2.1[) is dictated by a piecewise smooth referential model surface S. 
Imitating the classic approach of Kirchhoff, we specify the admissible kinemat- 
ics (Section I2.2j) by choosing for the displacement field in the tubular region 
Q(S,e) a representation parameterized by a few fields defined over S. In the 
deformations we envisage, (i) thickness may change, if the applied loads require 
and the boundary conditions permit; (ii) material fibers orthogonal to S must 
remain orthogonal to the deformed shape of S itself, an internal constraint we 
refer to as unshear •ability. Both the balance and the constitutive equations of 
our shell theory (Sections 12.31 and 12.41 resp.) are inherently consistent with the 
corresponding equations of three-dimensional linear elasticity: 

- the balance equations follow from a two-dimensional Principle of Virtual Pow- 
ers that is a direct consequence of stating the corresponding three-dimensional 
Principle for all virtual velocity fields in the linear space to which the admissible 
displacement belong; they are expressed in terms of a pair of two-dimensional 
stress measures that are defined as weighted thickness averages of the three- 
dimensional stress field in Q{S,e) (Section [23)1 : 

- the constitutive equations are arrived at when the three-dimensional constitu- 



2 




Figure 2: Shell axis is chosen parallel to roll- up axis. 

tive equations for unshearable orthotropic materials are inserted in the defini- 
tions of the two-dimensional stress measures. 

The remaining part of the paper is dedicated to cylindrical shells. We begin 
with shells whose thickness can change. In sections from 3.1 to 3.4, we parallel 
and specify the developments of Section 2 as to, respectively, geometry, kine- 
matics, balance laws, and constitutive equations. Then, we confine attention to 
axisymmetric equilibrium problems, and solve explicitly and exactly those of tor- 
sion and axial traction (Sections 4.1 and 4.2) - the cases for which experimental 
tests and numerical simulations seem to be especially easy to set up for CNTs 
- as well as the problems of pressure and rim flexure (Sections 4.3 and 4.4). Fi- 
nally, we lay down natural geometrical notions of thinness and slenderness, and 
we show how remarkably the analytical solutions derived in the previous section 
simplify for slender shells (Sections 5.1 and 5.2), and how effective contraction 
moduli and effective stiffnesses can be defined for a cylindrical shell, regarded 
as a traction or torsion probe (Sections 5.3 and 5.4). 

Next, in Section 6, we take up orthotropic shells whose thickness is consti- 
tutively immutable, a class that we designate by the names of Kirchhoff and 
Love by analogy with the corresponding classic plate theory. We adapt to the 
simpler case of Kirchhoff-Love cylindrical shells all the formulas derived in Sec- 
tions 4 and 5, both for whatever thinness and in the small thickness limit; in 
the latter case, we show how the four constitutive moduli characterizing the me- 
chanical response could be determined on the basis of simple real or computer 
experiments. 

In our final Section 7, we briefly recapitulate our main findings, and we 
indicate the directions of our future research, with special attention to the ap- 
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plication to CNTs of the concepts and methods developed in the present paper. 

2 General Theory 
2.1 Geometry 

In this opening section we recapitulate some well known notions, with the main 
purpose of introducing our notation and terminology. 

Following the approach to construct a shell theory proposed in [TO], we let 
S denote a compact, regular, orientable and oriented surface embedded in the 
three-dimensional Euclidean space £, and we let x denote its typical point and 
n(x) the value of its normal vector field at x, with |n(x)| = 1. We choose an 
origin o £ £, and denote by a; := x — o the position vector of x with respect to 
o. We assume that S admits a tubular e— neighborhood Q(S, e) (see Section 2.2 
of [5]) and a global parametrization 

R 2 D S 3 (z\z 2 ) i-> xiz 1 , z 2 ) £ S c £ 

(here S is an open set). A point p £ Q(S,e) has position vector 

p := p — o — x — o + (n(x), x £ S, ( £ I := (— e, +e), 

with respect to o; |£| is the distance of p from x, the point where the straight 
line through p perpendicular to S intersects S itself. The mapping 

(z\z 2 , C) ^ p(z\ z 2 , C) := x{z\ z 2 , C) + (n(x(z\z 2 )) 

is a global parametrization of G(S, e), with (z , z 2 , Q the triplet of normal curvi- 
linear coordinates of p (Figure[3]). We term the region G(S, e) of £ a shell-shaped 




Figure 3: Geometrical equipment of a typical shell-shaped region. 

region, of model surface S and constant thickness 2e. The thickness of Q(S,e) 
can be visualized as the length, whatever x £ S one picks, of the material 
fiber F(x) through x perpendicular to S; clearly, T(x) :— {p £ G(S,e) \ p = 
x + Qn(x). 
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The chosen parametrization induces a system of coordinate curves on 5, 
described by the mappings z a i-> x(z 1 ,z 2 ) (a = 1,2). The tangent space 
T x ■= span{e Q (x)} to S at x = (z 1 , z 2 ) is spanned by the tangent vectors to 
the coordinate curves at that point: 

e a (z , z ) : — x a (z , z ). 

On taking the normal field to be n(x) = vers(ei(x) x e2(x)), and on setting 

q(x) :— ei(x) x e 2 (x) ■ n(x) > 0, 

the covariant and contravariant bases at x are, respectively, {ei(x), e 2 (x), n(x)} 
and {e 1 (cc), e 2 (x), n(x)}, where 

?(x)e Q (a;) := (— l) a n(x) x e Q+ i(a;) (modulo 2, a not summed). 

For the rest of this subsection we leave the indication of the typical point x e S 
tacit. Accordingly, we write 

P := &i <g> e l = e* ® e, 

for the metric tensor and 

S P := e Q ® e Q = e a <g> e Q 

for the surface metric tensor. 

A vector field w defined over 5 can be represented both in the covariant basis 
and in the controvariant basis: 

v = Vie 1 = v 3 ej, with Vi := v ■ a, v J := v ■ e J , = 1, 2, 3). 

Analogously, a second-order tensor field T can be represented as 

T = T^e, <g> e, = T„e 4 <g> e J = T]e, <g> = T/V <g> e 7 , 



in terms of its covariant, contravariant, or mixed components T y = T • e % ® e- 7 , 
T y = T • e, <g> e j; or T] = T • e l <g> ej and T/ = T • e 4 <g> e-? (i, j = 1, 2, 3). 

Remark. The physical dimensions of covariant and contravariant basis vectors, 
and hence of the corresponding components of vectors and tensors, may differ. 
To circumvent this difficulty is easy, whenever it so happens that 

e x ■ e 2 = 0. 

Simply, one introduces the so called physical basis at x, that is to say, the 
orthonormal basis 

e e a 

{e<i>(x), e<2>(a;), n(x)}, e< a > := — — = — . 
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We shall be making use of physical bases and components (e.g., v<h> := v ■ 
e<h>) in Section [3l where we deal with shells whose model surface is a right 
circular cylinder. 

At a point of G(S, e), the covariant and contravariant basis vectors are, 
respectively, 

Qot ■■= P,a = e a + (n a , g3--=P : 3 = n; (1) 

and 

sg a := (— l) a n x g a+ \ (modulo 2, a not summed), g 3 = n, 

where 

C := Si x g 2 ■ n. 
The metric tensor G is defined to be: 

G := g l ® gi = g t ® g\ 

with 

S G := g a ® g a 

its surface part. 

The covariant bases {e^} at x G 5 and {^} at a; + C n ( x ) = p£ 5(5, e) are 
related by the s/iz/t tensor — briefly, the shifter — A: 

A{x, C) = 9i ® e i Ae k = g k ; 

it can be shown |10) that the ratio of the volume measures at p and x is equal 
to the determinant of A: 

dvol(p) = a(x, £) dvol(x), a := det A. (2) 

The shifter 

B(x, C) := g { ® e t = S B +n ® n, S B := # Q ® e Q , 

maps the controvariant basis at x into the controvariant basis at p. We have 
that 

B T A = P, AB T = G. 

Note that 

A(x,C) = s A(x,() + n(x)®n(x), A{x, () g a ® e Q = s P(x) - C W(a;), (3) 
where 

:= - s Vn = -n, Q (g)e Q 

is the curvature tensor of the oriented surface S. We have here denoted by 
s V the operation of taking the surface gradient of a smooth vector field v over 
S: S V« = v, a ®e a . Likewise, we denote by V the gradient of a vector field v 
defined over G(S, e): Vd = v,i®g l . Two divergence operators are associated 
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with the gradient operators, Div v :— V« • G and s Div v := S V« • S P, where the 
field v is defined, respectively, over Q(S,e) and over S. The surface divergence 
of a tensor field T over 5 is defined as follows: s Div( T T v) =: s DivT ■ v, for all 
constant vectors v. 

Let a the fiber- wise constant extension to G{S, e) of a vector field a defined 
over S. When taking the gradient of a, we have that 

Va = a, a ® ff Q = a a ® ( S B e Q ) = ( s Va) S B T . 

In the following, we will not make any notational distinction between a field 
defined over S and its fiber-wise constant extension to Q{S,e): e.g., we shall 
write: 

Va = ( s Va) s B T . (4) 

2.2 Kinematics 

The three-dimensional strain measure we use is the standard symmetrized gra- 
dient of the displacement field: 

E(u) = sym Vm := - (Vu + Vu T ). 

It follows from this definition that the covariant components of E are: 

2E ij = 2£ • ft <g> g 3 = (u k ® g k + g k ® • ( ft ® g 3 ) 

= {u ■ gt),j +(u ■ gj),i -u ■ (gi,j +gj,i), 

whence 

Eij = - (uij + u jti - u ■ (g hj +gj,i )J . (5) 

We restrict attention to shell-shaped bodies G(S, e) whose admissible defor- 
mations may induce thickness changes but must keep the material fibers or- 
thogonal to the model surface, in the sense that, whatever x £ S, the deformed 
material fiber u(T(x)) must be found orthogonal to the deformed model surface 
u(S). This pointwise internal constraint — a restriction on admissible displace- 
ment fields that, as anticipated, we refer to as unshearability — can be expressed 
in terms of the linear strain measure E in the following form: 

E(u)n ■ v = for all v such that v ■ n = 0, (6) 

or rather, equivalently, as 

E(u)n-g a =0 in G(S, e); (7) 
with the use of |TJ and (|147|> . © becomes: 

O • 9a), 3 +(« • n), a -2 n • n, a = 0, (8) 
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with the same quantification. We look for solutions of this system of two PDEs 
having the following form: 



u(x, C) = u (x) + C « (a;); (9) 



note, in particular, that 



u (x) ■ n(x) = E(u(x)) ■ n{x) ® n(x) 

is the stretch of the material fiber F{x), uniform all along it. Substituting ^ 
into ([5]), and exploiting the quantification with respect to the variable £, we find 

(o) (i) 

the following restrictions on the choice of the fields u , u over S: 

u ■ e a = — u , a n, [ u ■ n) , a = 0. (10) 

The second of (fit))) tells us that, for ([7]) to have a solution of type the fiber 

stretch must in fact be uniform all over Q(S,e); moreover, the tangential part 

(i) . (i) 

of u must be expressed as follows in terms of u and the curvature tensor of iS: 

« \ Sr7 /(°) \ „ r (o) 

u — [ u ■ n) n = — V ( it • n) — W u . 

With the use of ([3]), we conclude that we should pick: 

u s (x, C) := s A(x, C)a(x) + w(x)n{x) + C( - s V«'(i) + in(x)) , (11) 

where we have set: 

(0) ,(0) (0) (1) 

a := u — ( u ■ n)n 1 w := u ■ n, 7 := u ■ n. (12) 

Note that the displacement field (fTTj) is parameterized by two fields over 5, the 
tangential vector field a and the scalar field w, and by the constant 70 Note 
also that, by taking 7 = 0, we recover the kinematic Ansatz typical of the linear 
theory of Kirchhoff-Love shells [10] : interestingly, the displacement field: 

url(x, C) = s A(x, C)a(x) + w{x)n(x) - ( s Vw(i), a(x) ■ n(x) = 0, (13) 

can be shown to be the general solution of the system of three PDEs that 
expresses the Kirchhoff-Love constraint, namely, 

E(u)n = in g(S, e) (14) 

(cf. ([?]))• Remark. That the fiber stretch 7 defined by ([TH]) ,, be constant at any 
point of G(S, e) to satisfy the second of (|10[) is a counterintuitive consequence 
of expressing the unshearability constraint in terms of the linear measure of 



2 A vector field a (a tensor field T) defined over S is termed tangential if it so happens 
that a(x) ■ n(x) = (Tn = 0). 
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deformation E(u). In fact, on adopting one of the standard exact deformation 
measures, namely, 

D :=^((V/) T (V/)-l), where f(p) := p + U (p), 
one finds that 

D = E + ^H T H, H:=\7u. 

Thus, the exact conterpart of ([7]) is: 

= Dn ■ g a = En ■ g a + i Hn ■ Hg a , 

where the term quadratic in |_ff| does not vanish in general. In particular, for 
u of the form ([9]) , we find that 

„ „ (i) (o) (i) (i) 
Hn ■ Hg a = u ■ u , a +Q u ■ u , Q ; 

hence, the second of (|TUj) is replaced by 

(u ■n + |«f), a = 0, 

and intuition rings no bell. 



2.3 Balance Assumptions 

Given a material body occupying a three-dimensional open and bounded region 
Q and a velocity field v over f2, the internal power expenditure over a part II of 
fi associated with v is: 

[ S-Vv 
Jn 

where S denotes the stress field in ft; and, the external virtual power expenditure 
over II is: 

d - v + c - v, 
Jn Jdn 

where (d Q1 c Q ) denote, respectively, the distance force for unit volume and the 
contact force per unit area exerted on II by its own complement with respect 
to Q and by the environment of the latter. These representations of power 
expenditures are those typical in the theory of the so-called simple material 
bodies. In that theory, a part is customarily a subset of non-null volume of 
f2 (which makes it for a part collection deemed sufficiently rich), and a virtual 
velocity field is a smooth vector field whose support is a part; the Principle of 
Virtual Powers is the stipulation that 

S Wv = d ■ v + c v (15) 
Jn Jn Jan 



9 



for all virtual velocities in a collection modeled after a chosen collection of 
admissible motions, an alternative stipulation being 

IS Vv = d - v + c -v (16) 
Jn Jn Jan 

for all parts II of il and for all velocity fields v consistent with the admissible 
motions. 

The virtual-power principle is interpreted as a balance statement for the 
internal and external fields S and d Q , c a . As exemplified in [TT] for beams and 
plates, we use here below a restricted version of the three-dimensional statement 
(jTHJl to deduce a two-dimensional Principle of Virtual Powers and the associated 
balance laws for shells in terms of two-dimensional stress measures and applied 
loads. The first restriction we introduce has to do with the special shape of the 
three-dimensional bodies we consider: they must be shell-shaped in the sense 
of Section 12.11 The second restriction stems from the special class of admissible 
body parts we choose: they all must have the same thickness of the shell-like 
body. Thirdly and lastly, we pick a special class of virtual velocities, consistent 
with the representations of admissible displacements discussed in the previous 
section. 

2.3.1 Internal power expenditure. Stress measures 

The first two restrictions we listed imply that a typical part of Q = G(S, e) can 
be identified with the Cartesian product II = V x / of V, an open subset of S, 
and the interval /. As to virtual velocities, we choose them of the form ([5]): 

v(x,0 = ( v(x) + C { v(x), (17) 

with V the intersection of the supports of the vector fields v and v '0 On 
recalling (U]), we have that: 

Vt>=(V$ ) + C , V ( #)"B T + ( i ) ®n. 

Hence, 

S • Vd = {S S B) ■ B V ( v + ((S S B) ■ s W { v + Sn ■ %\ 

so that, having recourse to Fubini theorem and recalling ([2]), the internal power 
expenditure reads: 

/ S-Vv= [ ( 8 F- s Vu ) + 8 M- s V ( i ) +/ (3) • { v). 
Jn Jv 

3 Choosing instead virtual velocities of the less general, constrained form would pre- 
clude the appearance of reactive terms in the balance equations following from (1 16 I t . We shall 
return on this important issue after completion of the generating procedure we are now about 
to start, in Remark 12.51 at the end of Section 12.51 
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Here we have made use of the following definitions: 
s F(x) :=J i a(xX)S(xX) s B(xX)d(;=(J i a(xX)S(xX)9 a (xX)d()®e a , 

s M(x) := J a(x,CKS(x,() s B(xX)d( = ( / a(x,OCS(x,C)g a (x,QdC) ® *x, 

(18) 

and 

fW(x) >.= Ja(x,OS(x,Qn(x)dC, (19) 

for the two-dimensional stress measures that we call, respectively, force tensor, 
moment tensor, and shear vector. 

It is important to observe that the force and moment tensors are tangential. 
Now, for a tangential tensor field T over S, it follows from the general identity 

T ■ s Vv = s Biv(T T v) - s DivT • v 

and the standard divergence theorem that 



T- s Vv= / Tmv- s DivT-v, 

V JdV JV 

where, if t is the tangential vector to the curve dV in the tangent plane to the 
surface S, m = t x n is the normal to dV . Thus, 

/ SVv= [ (- s D\v s F- ( v+(- s mv s M+f^)- { v)+ [ ( s Fm- ( v+ s Mm-v 
Jn Jv ^ ' Jdv ^ 

2.3.2 External power expenditure. Applied loads 

As to power expenditure of the applied forces, we find: 

J d a - v = J {{J ad o) ■ ( $ + ( J a(d ) ■ t>^ 
for the distance force, and 

/ c a v= f(a + c+ + a~c~) ■ v + e(a + c+ - a~c~) ■ V 

JdU JV V 

for the contact force0 where we have set: 

a ± (x) := a(x, ±e), c^(x) :— c a (x,±e), x € V. 

We are now in a position to define the load fields over S induced by the three- 
dimensional fields (d , c ). These are: 



(20) 



4 Here we have made use of the fact that dll = {V X {±e}} U {dV X /}. 
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• the distance force and distance couple per unit area 

q a (x) := / a(x,()d (x,()d( + a + (x)c+(x) + a~(x)c~(x), 

Jl f (21) 
r a (x) := J a(x,Q(d (xX)d( + e(a + (x)c+(x)-a (x)c a (x)); 

• the contact force and contact couple per unit length 

l (x) := J c Q (x, ()d(, m (x) := J^c (x,()d(. 

All in all, the external virtual power expenditure takes the following form: 

d - v + c ■ v = / (q ■ v + r a ■ v ) + / (l ■ v + m ■ v ). 
Jn Jon Jv Jdv 

2.3.3 Principle of Virtual Powers. Field equations. Boundary equa- 
tions 



The two-dimensional Principle of Virtual Powers we arrive at is 
0= / ((- S DW S F -q a ) ■ { v + (- s Div s M+/ 

r o i (22) 

+ J ({ s Fm-l )- { v +( s Mm-m„). ( i ) ), 



- ro) ■ (1 v 



(o) (i) 

a statement to hold for every pair of vector fields v and v on S and for every 
part V of S. Under the standard blanket assumptions of smoothness, and with 
the use of a standard localization lemma, (|22l) yields the (two-dimensional) field 
equations to hold at any interior point of S: 

s Viv s F+q o = 0, 
s Div s M-/ (3) + r a = 0. 

Granted (1231). what remains of 



/ f ( S F m-l )- { v + ( S M m - m ) ■ ( v) = (24) 
Jdv v ' 

(0) (1) 

for all admissible variations v and v . Localization of (|24j) yields different 



results according to where it is performed. At an interior point of 5, where it 
can be combined with arbitrariness in the choice of V, we have that: 

s Fm = l , s Mm = m Q , (25) 

two relations that parallel the classic relation between stress tensor and contact- 
force vector for three-dimensional Cauchy bodies^ At a point of dS where a 

5 See |15) for a discussion of this issue that covers second-gradient materials as well. 
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Dirichlct boundary condition prevails - that is to say, where one or more com- 
ponents of the boundary trace of the displacement field are prescribed - the cor- 
responding components of the admissible variations must vanish; accordingly, 
the complementing components of both vectors ( s Fm — l a ) and ( s Mm — m ) 
must also vanish (yielding boundary equations of Neumann type), because lo- 
calization of (1241) leads to: 



( s Fm - l ) ■ ( v = and ( s Mm - m ) ■ ( v = 0, (26) 
(o) (i) 

for all admissible choices of v and v . 

Remark. The shear vector f( 3 > , defined by (fH))) and entering the last of equations 
(|3"2"|) , is one of the force vectors : 

f®(x) := Ja(x,QS(x,Qg\x,Qdt; (27) 

the moment vectors are: 

m«(z):= [ a(x,()(S(x,()g l (x,()d(. (28) 



With these definitions, we may further set: 

F := / w ® ej = a F +/ (3) <g) n, M := m w ® e 4 , 

with 

s ^ = ® ec, S M = ® e Q , (29) 

and 

f( a )= s Fe a , m (a) = s Me a . (30) 
In terms of force and moment vectors, equations ((23)) become: 



/ (tt) , Q +7 Q V (a) + 9o = 0, 



m 



i8- 

'" ) -+l^m {a) -f i3) +r o = 0. 



(31) 



We call membrane forces the components := F ■ e a ® e' 3 — J'' 9 ) • e Q = 
S F -e a (£) e@ - respectively, normal m. f. for a = f3 and shear m.f. for a. ^ f3; 
and we call M aa :~ M ■ e a e Q = n»W • e Q = s M-e Q <g) e Q (a = 1,2) 
the bending moments and Af Q/3 := M • e a (g) = m (a) • e* 3 = S M ■e° t ® e 13 , 
(a, (3 = 1,2, a ^ ft) the twisting moments. Finally, we call F 3a — ■ e a the 
transverse shears, F 33 — f^ 3 — f^-n the thickness shear, and M 3a = m^-e a 
the thickness moments. Component-wise, the general field equations (|23j) can 
be written in the following form: 

F Sa \ a -W^F 3a + q 5 o =0 ((5 = 1,2), 
F 3a \ a + Wf, a F^+q 3 o =0; 

M Sa \ a - F 3S + r 5 a =0 (<5 = 1,2), 
M 3q | q + Wp a M t3a - F 33 + r 3 Q = 0, 



13 



where a vertical bar denotes covariant differentiation (given a tensor field T 
over S, T Sa \ a := Tf° + lp a T ?OL + 7^T Sa , with T 5a := T ■ e 5 <g> e a ) and where 
ji a :— e 5 ■ ep, a are the surface Christoffel symbols. 

Remark. A shell-shaped body Q{S,e) is in a membrane regime if it so happens 
that the following conditions: 

/(<*) . n = 0, / (3) = 0, m< Q > = (33) 

hold identically in 0(<S, e). The first two of fl33]) imply that F 3a = 0, the fourth 
that M la = 0. Consequently, equations (1321 reduce to: 

w^f** + = 0, 

together with the following compatibility condition on the data: 

r Q = 0. 



2.4 Constitutive Assumptions 

The penultimate step in our construction of a shell theory consists in specifying 
how the two-dimensional stress measures fTS]) depend on the admissible defor- 
mations pip . We wish to come up with the simplest theory accommodating 
both an orthotropic response and the unshearability constraint discussed in Sec- 
tion 12.21 To this end, we confine attention to cases when the material response 
is uniform all over G(S, e). Moreover, all along each fixed material fiber !F{x) we 
take the orthotropy plane orthogonal to that fiber, so that the orthotropy and 
the tangent planes coincide at x £ S; in particular, whenever use is made of or- 
thogonal curvilinear coordinates, as is the case for the cylindrical shells treated 
in Section [31 we take the orthotropy axes tangent to both the coordinate lines 
z a = const. 



2.4.1 Orthotropic elasticity tensors 

Let C denote the fourth-order tensor of elasticity, a linear transformation of the 
space of symmetric tensors, that specifies the stress response to deformations 
in the parent three-dimensional theory we are going to select. An orthogonal 
tensor Q is a symmetry transformation for the linearly elastic material described 
by C if it so happens that 

QC[E] Q T = C[QEQ T ], for all symmetric tensors E, 

i.e., that 

QC = CQ, 

where the fourth-order tensor Q is such that Q[A] := QAQ T for every A in 
the space the second order tensors Lin; the symmetry group of C is the set 

^ C :={QeOrth|QC = CQ}. 
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Given a subgroup Sf of the orthogonal group (or, for what it matters here, of 
the group of all rotations), one seeks a representation formula for all elasticity 
tensors C such that £fc 3 ^, i.e., for all elasticity tensors sharing a given sym- 
metry group. A linearly elastic material is called orthotropic when its stress 
response is insensitive to a rotation of 7r about a given axis c, i.e., when the 
symmetry group of its elasticity tensor includes that rotation. We give here 
below a general representation formula for the elasticity tensors in question. 

Let {ci (i = 1,2,3)} be an orthonormal basis of vectors. Consider the 
following orthonormal basis for the linear space Sym of all symmetric tensors: 

V a = -^=(c a ® c 3 + c 3 <g) c a ) (a = 1, 2), V 3 = c 3 ® c 3 , 

^ ! (35) 

W„ = c Q ® c a (a not summed), W3 = —= (c% ® c 2 + c 2 ® ci). 

With the use of this basis, any orthotropic elasticity tensor can be written in 
the following form: 

C = dm Wi ® Wi + C2222 W2 ® W2 + C3333 V 3 ® V 3 + 
+ C1212 W 3 ® W3 + C3131 Vi ® Vi + C2323 V 2 <E)V 2 + 
+ Cii 22 ( Wi ® W 2 + W 2 ® Wi) + Ch 33 ( Wi ® F 3 + F 3 ® Wi)+ 



+ C 2 233( W 2 ® F 3 + F 3 ® W; 



2 1 



(cf. [B]); the orthotropic material class is then parameterized by the 9 elastic 
moduli Gnu, ■ • ■ ,C 2233 . 



2.4.2 The elasticity tensor of unshearable orthotropic materials 

As a direct consequence of the fact that the shell model we are after incorporates 
a kinematic constraint, the appropriate three-dimensional response is captured 
by an elasticity tensor somewhat simpler than (|3^|) . To see why, and to derive 
such an elasticity tensor, we apply to our present case a general representation 
result in constrained linear elasticity [TH [11] . 

We make the shell geometry agree with the geometry intrinsic to the material 
response, in the sense that, at any fixed point x £ S, we identify c 3 with 
n(x). With this identification, the internal constraint (JB]) can be read as the 
requirement that all admissible strains be orthogonal to the following subspace 
of Sym: 

2/ 1 - := span(F Q , a = 1,2). (37) 

Accordingly, the space Sym is split into the direct sum of two orthogonal sub- 
spaces: 

Sym = ^®^ ± , (38) 
and the stress is split into reactive and active parts: 

S = S {R) + S {A \ (39) 
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with 

5 (A) =C[JB], C:@^3>; 

here the coefficients ^„ are constitutively unspecified, and the constraint space 
implicitly defined by (|37|) and (j38|) can be identified as 

^ = span(V 3 ; W u i = 1,2,3). (41) 

On applying a general result proved in [14] , a representation for the desired 
elasticity tensor C can be deduced from the one given for C in (|36|) : 

C = P@C|@, P@ :=I- V Q ® V a . 

where P® denotes the orthogonal projector of Sym on 3>. One finds the 7- 
parameter representation 

C = dm W l ®W 1 + C2222 W2 ® W2 + C3333 V3 ® V3 + C1212 Ws ® W3+ 
+ C U22 ( Wi ® W 2 + W2 ® Wi) + Cii 33 ( Wi ® V3 + V3 ® Wi)+ 

+ C 2 233( ® F 3 + V 3 ® W2). 

(42) 

As is customary, we assume that C is positive-definite, i.e., that 

E ■ C[E] > for aH E G 9 \ {0}. (43) 

Remark. In case the internal constraint (j6]) is reinforced a la Kirchhoff-Love as 
specified by (fl"4"|) . so as to exclude thickness changes, the above procedure yields 
the elasticity tensor 

C = dm Wi® W1+C2222 W 2 ® W2+C1212 W3® W 3 +Cn22( Wi® W 2 + W 2 ® Wi), 

(44) 

a linear transformation of the constraint space 9 = span(Wi, i = 1,2,3) into 
itself. 

2.4.3 The technical moduli 

In technical applications of classic isotropic elasticity, the two Lame constants 
are replaced by the technical moduli E and v of Young and Poisson, plus the 
shear modulus G, under the condition that E = 2(1 + v)G. Likewise, we here 
replace the seven Lame-like constants dm, . . . , C2233 in (|42|) by an equivalent 
list of ten technical moduli - three of them being Young-like, six Poisson-like, 
and one shear-like - that must satisfy three independent algebraic conditions. 
The technical moduli in question are precisely those that enter the following 
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representation of the compliance tensor 



-l. 



cr 1 = — Wi ® Wi + ^- w 2 ® w 2 + — v 3 ® v 3 + — ^— w 3 ® w 3 + 

- ^( Wi <g> W 2 + W 2 (8 Wi) - ^( Wi «) F 3 + V 3 (8 W x )+ (45) 

- v -^(w 2 ® v 3 + v 3 ® w 2 ) 

E 2 

(as is well known, the positivity assumption (|43p guarantees invcrtibility of C). 

To confirm that the three moduli Ei,E 2 and E 3 are Young-like, imagine 
to induce a state of uniaxial traction in the direction, say, C\ in a specimen 
made of the material under examination, so that the stress is S = a W\ and 
the corresponding strain is 

E = aC- 1 [W 1 ] = |-(Wi - v 12 W 2 - u 13 V 3 ). 
Then, we find that 

the ratio between the axial stress and the corresponding axial deformation. 
Moreover, we find that 

E -W 2 E-V 3 
Vl2 = -E-Wl' ^ = -E-W^ (4?) 

the Poisson-like negative ratios of the transverse deformations in the directions 
c 2 and c 3 and the axial deformation in the direction C\ . The defining formulae 
for E 2 , i>2i, v 23 and E 3 , v 3 \, v 32 are completely analogous to (|46l) and (|47l) . Due 
to the built-in symmetries of C, it turns out that 

Ei _ v\ 2 Ei _ v\ 3 E 2 _ v 23 
E 2 v 2 \ E 3 v 3 \ E 3 v 32 



Finally, the formula 



2G= S - W * 



E-W 3 

establishes the nature of G as a shear modulus. 

It is straightforward to see that G — C\ 2 \ 2 . The other technical moduli can 
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be written as follows in terms of the components of 1 
£1 



C1111C2222C3333 — 


C3333C 1122 - C 2222 Cf 133 


-c im « 


~<2 

-2233 


h2Cn22<I 


-1133C2233 




C2222C3333 — c 


i2 
'2233 








C1111C2222C3333 — 


C3333Ci 122 - C2222Ci 133 


-Ci U1 « 


~>2 

-2233 


V 2Cll22« 


-1133C2233 




C1111C3333 — C 


'1133 








C1111C2222C3333 — 


C3333Ci 122 - C 2 222Ci 133 


-Cunt 


~>2 

-2233 


K2Cii2 2 d 


-1133C2233 



E2 = 

E 3 = 



C1111C2222 — C? 



1122 



3333^1122 — C1133C2233 C3333C1122 — ^1133^2233 



C2222C3333 — C 2233 ' C1111C3333 — Ci 133 



^2222Cn33 — Cii22C 22 33 C2222C1133 — Cn 22 C 



2233 



C2222C3333 — C 2233 ' CiinC 2222 — cf 122 

_ CimC 22 33 — C1122C1133 _ CiinC 22 33 — Cn 22 Cii33 

rn in ip2 ' ^32 IT (T (T"2 

^1111^3333 — ^1133 ^1111^2222 _ ^1122 

(49) 

To sum up, the 7-parameter representation (|42|) can be re-written in the 
form: 

C = A _1 ^£'i(l - l>23 ^32) Wi <g> Wi+S 2 (l-^i3^3l)W" 2 (8) W 2 

+ E 3 (l - u 12 V21) V 3 ® V3 + 2A G W3 W 3 

+ S x (z/ 21 + i> 23 1/31 )( Wi (g) W 2 + W 2 Wi) (50) 
+ Etfai + V21 vz%){ Wi ® V 3 + V 3 ® Wi) 

+ E 2 {v 3 2 + V 31 Vi2)(W 2 ® V3 + V3 ® W 2 )J, 

where 

A := 1 - z/ 12 // 21 - z/ i3 i/ 31 - ^23 z/ 32 - 2vi2 v 23 v 3 i- (51) 

2.5 Force and Moment Vectors and Tensors 

The final step in the assemblage of our shell theory - the posing of initial- and 
boundary-value problems - demands that the balance equations (|32[) are writ- 
ten in terms of the parameters involved in the general representation (|11[) for an 
admissible displacement field. All we have to do to construct the corresponding 
parametric representations for those components of the force and moment ten- 
sors that enter the balance equations is to insert in the definitions (|27|) and (I28|) 
the general parametric representation for the stress field in Q(S,e), and then 
make use of definitions (|29[) . Now, the required stress representation is obtained 
by combining ([39 ]) -(|30 ]l - with c 3 = n(x) - and (g2j - with E = E(u s ), and 
v,g given by (jlll) ; one finds: 

S(x,C) = s<- R) (xX)®n(x) + n(x)® s^ R \x,Q +C[E(u s (x,Q)}, (52) 
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where the restriction to the fiber F(x) of the reactive field s^ R ' is a vector field 
perpendicular to n(a;)[f| Consequently, the force and moment vectors have both 
reactive and active parts, namely, 

/ (Q) (z) = (J a(x, 0(s (R) (x, • 9 a {^ C))dC) n(x), 

f (3) (x) = J a(xX)s {R) (x,C)d<:, = J a(x,OC[E(u s (x,OW(x,C)dC, 

(53) 

and 

m^(x) = (jf a(z, 0C(* (fl) (*, • 9 a (x, C)R) n(x), 

m^(x) = J a(xXKsW(xX)d(, m®(x) = J a(^ C)CC[£(«s(x, C))]<7 WR- 

(54) 

It is not difficult to check that 

f (i) (x)-f^(x)=0, m (l) (x) -m {l) (x) =0. (55) 

It is also easy to see, in the light of (Pit?)) , that the force and moment tensors S F 
and S M have active and reactive parts as well. Thus, the balance equations one 
arrives at are not pure, in the sense that, in addition to the parameter fields, 
they also include reactive terms. In fact, the transverse shears are reactive, and 
equations ([3"2"j) 3 relate them to the active bending moments: 

F 35 = M Sa \ a +r s R (56) 

The thickness moments are also reactive, and (|32[) 5 relates their divergence to 
the active bending moments and the active thickness shear: 

M 3a \ a = -Wp a M^ a + F 33 -r 3 . (57) 

These observations suggest a sequential strategy to solve a shell problem 
within our present theory, where the unknowns are (the fields a,w and the 
constant 7 that parameterize) the displacement field us and the reaction force 
and moment fields: firstly, by the use of the projection operator defined 
in the Subsection 12.4.21 one derives a set of reaction-free consequences of the 
balance equations; secondly, one solves such 'purified' system of equations for 
us', thirdly, one returns to the full balance equations, where the active terms 
can now be computed explicitly, and solves them for the reactive fields H 

6 It follows from the first of (gOj and the first two of J35]) that = (g) c 3 + c 3 (g) , 

With *(«)= Jy^Cc. 

7 In fact, F 3S = /< 3 ) ■ e 5 = /< 3 > ■ e s ; and, M Sa = mW ■ e a = mW ■ e a . 
8 The annihilation procedure of the reactive terms occurring in the three-dimensional bal- 
ance equations of constrained linear elasticity is discussed in Section 17.2. 
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Remark. As anticipated in footnote 2, reactive forces and moments are found in 
the balance equations (|32l) , because the class of variations (fT7|) we used to derive 
those equations from the Principle of Virtual Powers (|16[) is visibly larger than 
the class (|TTj) of admissible displacements (this will not be the case in the next 
section). Although we here do not pursue this issue any further, we recall that 
having at one's disposal the reactive dynamical descriptors associated with the 
kinematical Ansatz adopted to construct a lower-dimensional structure theory 
can be proved beneficial to improve the pointwise approximation of the relative 
three-dimensional stress field [SI El [El E2 • 
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3 Cylindrical Shells: Generalities 

3.1 Geometry 

We now restrict our attention to shells whose model surface S is a portion of 
a right circular cylinder, that we parameterize as usual by means of cylindrical 
coordinates: 

z 1 — x\, z 2 — d 

(Figure 2]). For p Q the radius of the directrix of S and n the normal to S, and 



Figure 4: Geometrical equipment of the model surface of a cylindrical shell. 

for e < p a the thickness of the shell-shaped region Q(S, e), the position vectors 
with respect to the origin o of two typical points p g Q(S,e) and x G S are, 
respectively, 



x — o = x(xi,-d) = x\C2 + p n{d), n($) = sin^ci + cos$C3. 

The relative covariant and contravariant bases can both be represented in terms 
of the orthonormal basis {ci, n' (■#), n(i?)}, that serves as physical basis (recall 
Remark l2~Tj) : 





PoJ 



and 



9i(p) = ei{x), ei(x) = a, 




(58) 



53 0) = e 3 (x) 



n{§). 



and 



g 1 (p) = e 1 (x), e 1 (x) = cj, 




(59) 



g 3 (p) = e 3 (x) 



n(iJ). 
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The surface shifter defined by ([3]) 2 turns out to be: 

S A($, C) = ci ® a + (l + j-^j n'(<&) ® n'(0), 

whence 

a(() = 1 + ^ . 

Po 

Finally, the only nonnull Christoffcl symbols on S are: 

722 = ~Po> 723 = 732 = Po 1 - 

3.2 Kinematics 

The displacement field (fTTj) now reads: 

«c(xi,i?,C) = fci <g) Ci + fl+ -^J n'(tf) ® n'(i?)J a(a;i,i?)+ 

+ (w(a;i,<?) +C7)»W - cf«U(*i>0)d +p-S(si ) tf)n'(tf) 

(60) 

its representation in the physical basis is: 

uc{x\,d, Q = u<i>(xx, C)ci + u<2>(xi,i}, C)n'(i?) + u<3>(^i,i?, £)n(i?), 

(61) 

with components: 

u<l> = a<i> — (w : i , 

u<2> = | 1 + — | a<2> — —w 2 , (62) 

V Po/ Po 

u<3> = w + £7 . 

Since 

1 / ( \ _1 

V«c = « i ® = m,i (8 ci H 1H « 2 <8> «' + u 3 ® n, 



Po V Po 
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the nonnull physical components of the strain tensor are: 

E<n> = Vt*c • c\ ® c\ = u<i> j i = a<i> ! i — C w ,n > 
2E<12> = 2E<2i> — Vu c ■ (ci ® n' + n' ® Ci) = u<2> i + p Q ( 1 + — ] u<i> 2 

V Po/ 

1 + — ) a<2> i i — — iu 21 + — f 1 + — ) (a<i>.2 - C w ,i2) 7 

Po) Po Po V Po/ 

i?<22> = Vt*c ■ n' (8 rc' = p ( 1 H J (M<2>,2 + U<3>) 

\ Po / 



(63) 



a<2>. 2 U' 22 + w + C7 

PO V Po/ VV Po/ Po 

i?<33> = V«c ■ n ® n = u<3>,3 = 7 . 



3.3 Balance Assumptions 

We exploit the virtual-power procedure detailed in Section 12.31 with the differ- 
ence anticipated in Remark 12.51 the variations we now employ have the same 
structure as the admissible displacements (|6Tj) - (j62j) : hence, no reactive contri- 
butions are going to enter the field and boundary equations that the procedure 
delivers. Precisely, the variations in question have the following form: 

(°) / , 

v = v\Ci + v 2 n + v 3 n, 

(i) _ (64) 

v = -v S)1 ci + p a («2 - v 3 , 2 )n' + v 4 n, 

with the scalar fields Vi — Vi{x\,$) (i — 1,2,3) compactly supported in S and 
with W4 a constant. 

3.3.1 Field equations 

With a view toward deriving the field equations to which the general balances 
(|23|) reduce in the present situation, we firstly take w 4 = in (|64|) . In this 
instance, on setting: 

(o) fol 

/ := -( s Div s F+q ), m := -( s Div s M -/< 3 > + r ), 
the formulation (|22[) of the Principle of Virtual Powers reads: 



- f M (0) , (0) (1), f ( M , M _! (0) 

= J (f ■ v + m • v) = J [vx ( / • ci) + v 2 ( / + p a 1 m) • n + 



S°) /(0) \ ,(o) A . 

+ v 3 {f + (to- aja+Po {m-n), 2 )-n 
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Hence, three field equations must hold at each interior point of 5, namely, 

(0) 

/ • C! = 0, 

(/ + /^)-n' = 0, «*) 

(/ + (m ■ cij.i +p Q [m ■ n ), 2 ) ■ n = 0. 

To find the balance law that follows from testing the equality of internal and 
external powers by way of virtual velocity fields of the form: 

v = Q v , v = V4n, (66) 

with V4 an arbitrary real number, we turn to the formulation (fTS)) of the Prin- 
ciple: indeed, V4 must be taken constant over the whole of fi, because it is 
meant to be a variation of the constant thickness strain _E<33> = 7 in the rep- 
resentation (|62p of the admissible displacement^ Given that, for any v as in 
(EH), 

Vw = V4(n (g) n + — n' <£> n'), 

Po 



and that V4 can be chosen arbitrarily, (fl5j) and ([66)1 yield: 

/ / if ^ ' n<g> n+ P~ 5 ~ ^ n ^ d ^) dXlM = ( ' 67 ' ) 

In terms of physical components of the force and moment tensors, the system 
of equations (j65|) and (l67l) can be written as: 



F<n> 1 + p o 1 F<12>, 2 + <?o<l> = 0, 



(T<2i> + p a 1 Af<2i>) ! i + p a 1 (F<22> + p Q 1 A/<22>)^ 2 + g <2> + P a 1 r a <2> = 0, 

M<11>,H +/9o 1 (A/<12> + A/<21>) ! 12 + \M<22>.22 + 

Po 

- p^F<22> + q a <3> + r D <i>4 + p~ V Q <2>, 2 = 0, 

/ 1 A/<22> + F<33> — r Q <3>) dxidtf = 0; 

-z Jo 

(68) 

9 More generally, we point out that, whenever the admissible motions are such as to suggest 
the use of test fields of the form v = v w, with v a an arbitrary constant and w a given vector 
field, then l|15H is the appropriate formulation of the Principle of Virtual Powers, and the 
associated balance information is an integral relation of the form: 



S Vwi = d ■ w + / c 

Jn Jn Jan 
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the following version in terms of contravariant components may easy a compar- 
ison with ([32|) : 

F 1 \ 1 +F 12 l2 +q 1 o = 0, 

F 21 ;1 +F 22 ll +p- 1 (M 2 \ 1 +M 22 , 2 +r 2 ) + q 2 = 0, 

M u , n +M 22 , 22 + (M 21 + M 12 ), 12 - Po F 22 + q 3 a + r 1 ^ + r 2 ^ = 0, 

I (p M 22 +F 33 - rl) dxidtf = 0. 
-i Jo 

3.3.2 Boundary equations 

We start from the general weak statement (|2"1|) . repeated here for the reader's 
convenience: 

/ (( S F m-l )- { v +( s Mm- m ) ■ ( v) = 0, 

where we now insert variations of type (|64p . For simplicity, we restrict attention 
to parts of the model surface whose boundary consists of the union of two 
identical parts do, d\ of directrices at different axial abscissae and the relative 
two segments of generatrices go, gi, that is to say, with reference to Figure [SJ 
parts V such that 




Figure 5: A typical part of the model surface of a cylindrical shell, bounded by 
generatrices and directrices. 

dV = \Jd a Ug a (a = 0,1), 

a 

with 

da := {x eS| x = (Z 0) i?), d G (0 O) i?i)}, g a := {x G S \ x = (xwd a ), %i € (lo,h)}- 
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Over d\, where m — ci, the integral condition ([M|) takes the form: 



= / [{ S F d - l ) ■ (vici + v 2 n' + v 3 n)+ 

( S M ci - m ) ■ (-v 3 ,i ci + p~ x {v 2 - v 3 , 2 )ri + v±n)) 

J (i 7 '<ll> — 1 <1>) + V 2 (F<21> + p Q 1 M<21> — l a <2> — p a m <2>)- 



xi=h 



V 3 (F<31> - p o 1 M<21> ! 2 - / D <3> - p o 1 m a <2>^) + 



— u 3 ,i (M<u> — m <i>) + W4(Moi> — m <3>)^ 
+ (v 3 (M <2i> - p~ 1 m <2>fj 

over jji, where m = n'(d), becomes: 



(Ji,i?o) 



(69) 



0=/ (( s .Fn'- Z ) • (wici +u 2 n' + u 3 n)+ 

( S M n' - m ) • (-u 3 ,i ci + p" 1 ^ - «3,2 )«' + W4«) 
/ 



Wl(F<12> — Z <1>) + V 2 (F<22> + p o M<22> — l Q <2> — p Q m <2>) 
V 3 (F<32> — p~ M<12> t % — Z D <3> — p~ m <2>,x)+ 



— v 3 ^(M<22> — m D <i>) + Vi(M<32> — m <3>)j 
+ (v 3 (M<i2> - p~ 1 m <2>)j 



dx\- 



(Jo,t?i) 



(70) 

needless to say, relations completely similar to ([69]) and ([70)) hold, respectively, 
at do and go. 

. . . (1) (o) 

As mentioned in Section [51 if it 3 or anyone of the components of u and 

(0) 

V u 3 is assigned at a boundary point, then V4 or the corresponding component 

(0) (0) 

of v and V113 must vanish at that point; on the other hand, the arbitrary 
variation of each of the remaining parameters in (|64p induces at the same point 
a boundary equation of Neumann type: for example, inspection of ([69| shows 
(o) 

that, if u ■ C\ = a<i> is assigned over d a , then v\ must be taken there identically 
null, while the Neumann boundary equations 

i^<21> + p~ M <21> = l Q <2> + p^ > 1 m <2> etc. 

must hold; in other words, at a point of d a admissible boundary conditions must 
consist of a list of mutually exclusive assignments of the one or the other element 
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of the following five power-conjugate pairs: 

(F<u>,a<i>), (F <2i> + p^ 1 M<2i>, a<2>), 

(Foi> - 1 M<2i> i2 , w), (M<ii>,io,i), (Moi>,7). 

Likewise, from (|70[) we deduce that the boundary conditions at a point of g a 
should consist of mutually exclusive assignments of the pairs: 

(F<12>, <Z<1>), (F<22> + p~ 1 M<22>, a<2>, 

1 / \ ) 

(F<32> — p~ M<12> t i,W), (M<22>, W,2 ), (M <32>, 7). 

3.3.3 Inertial interactions. Evolution equations 

The time evolution of an unshearable orthotropic shell is ruled by the partial 
differential equations that follow from (|23|l when the inertial force is separated 
from the rest of the distance force per unit volume of G(S, e). To this effect, we 
set: 

d™ := -S ii and df := d - dfr, 

for, respectively, the inertial and noninertial distance forces (here 5 denotes 
the mass density per unit reference volume and a superposed dot signifies time 
differentiation). We then set: 



Qo n ■= - J <x(x, C) S o{x, Qu(x, C; t) d( and r^ n := - / a(x, 0( S o(x, C)«(a;, C; *) dC> 

(73) 

for the inertial force and the inertial couple per unit area of S. Finally, we 
return to definitions (pH]) -^ and set: 

q?~qo-q? and rf := r - r* 

for the relative noninertial loadings. In conclusion, the balance equations (|23l) 
take the evolutionary form: 

C = s Div s F+q l :\ 

ri n = s Div s M-/ (3) +rf. 
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In the case of cylindrical shells whose mass distribution is uniform, definitions 
CSJ yield: 

C<!> = -So (a<i> ~ \ £2 Po 1 * i 

- / / 1 e 2 \ le 2 



le 2 



r n <i> = -So^Po 1 a<i> - w t ij, 
- le 2 

r™<2> = -S -— (a<2> - ii) 2 ) i 
3p 



(75) 
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where (5 := (2e)5 is the uniform mass density per unit area of the model surface 
S. Hence, the evolution equations corresponding to the balance equations (|68[) 
are: 

So ^a<i> - 7^£ 2 Po 1, w,iJ = F<n>,i + Po 1 F<12> <2 + Oo'<l>, 

= (F<21> + / 9~ 1 M<21>) i i +p t 7 1 (F<22> +p t 7 1 M<22>) :2 + g ^<2> +p f 7 1 r^ i <2>, 

<5o + (7 + S<1>,1 - P W,11 + p~ 1 a<2>.2 - p^W^jj 

1 

= M<11> 11 + p" 1 (M<12> + M<21>) 12 + -rM<22> 22 ~ Po 1 F<22> + 

Po 

+ C<3> + C<1> 4 + p^rf <2>, 2 , 



So-Z W + Pol = Po M<22> + F<33> - r° l <3>. 

(76) 

Needless to say, these equations are to be equipped with a set of initial 
conditions for the unknown fields a<i>, a<2>, w, and 7, and for their time rates. 



3.4 Constitutive Assumptions 

The components of force and moment tensors that appear in (|68[) depend on 
the active part of the three-dimensional stress field. The latter has the following 
expression in terms of the constitutive law (l50l) and the strain tensor (|63[) : 

S(x,0 = C[E(uc(x,0)}; (77) 
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in components, this equation reads: 



EE E 

S<11> = —(1 - ^23^32) E<11> + "^-(^21 + ^23^3l) ^<22> + "^-(^31 + ^21^32) -E<33>, 

Eh E\ En 

S<22> = — (1 - V 13 V 31 ) E<22> + "^-(^21 + ^23^3l) E<11> + "^-(^32 + ^31^12j -E<33>, 

E^ E\ En 

S<33> — -^-(1 - V i 2 V 2l) E<33> + "^-(^31 + ^21^32) E<U> + "^-(^32 + ^31^12) E<22>, 
S<12> = S<21> = 2G E<12>, 

(78) 

where the constitutive modulus A is defined as in (f5Tj) . the components E<a/3> 
as in (|63p . On inserting these relations in the appropriate consequences of 
definitions (|2U)) . we find: 



F u = / ( l + -^]s<n> = 2 — ^ 

PoJ Po A 



(1 - 1/23^32) ( P a<i>,i - 3 e2w ,n 



+ {V2l + ^23^32)(a<2> j2 + w) + (2/3I + ^21^32)p Q 7 



F<22> = I S'<22> = 2— ^ 

Po A 



1 



1 , 1+ 

Po P. 



(1 - I/I3^3l) ^£K2>,2 - I 1 - l0g J __ 



!j- to ) + —(^21 + ^23^31) a<i>.i + (^32 + ^31^12)P Q 7 

J V 

(80) 



(79) 

(10,22 ~ Pol) + 



F<33> 



1 + 1) S<33> = 2 ±E<L 
PoJ Po A 



(1 - ^12^2l)p o 7 + 



^31 + ^21^32 
A 



1 2 

p Q a<i>,i - -e w t u 



F<2i> = J + > ; ^ 2,-r; 

and 



1 £< 



^32 ^ ^31^12 
M 



le 2 



(a<2> i2 + to) 
(81) 



1 + 77^ a<2>,i - --zW,2l H a<i>.2 

3 Po J 3Po 

(82) 



1 

Po 



M u = / I l + i-V5<xi>=~^ 



(1 - ^23^32)(a<l>,l - p o W ll) + 



+ (l/ 2 i + ^23^3l)(a<2>,2 - ^,22) + (^31 + ^21^32)7 



(83) 



C \ 4 e 3 

M<2i> = I I H ^5<2i> = - — G(a<2> i — to i 2 ), 

//V PoJ 3 p Q 



(84) 
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M<12> = / C5<12> 



a<2>.i — W 91 



1+ t 



1 



1 - log 



Pa 



1 - 

Po , 



10,12 



a<i>,2+ 



(85) 



M<22; 



where 



CS'<22> = 2e- 



(1 - ^13^31) 1 



■log 



Pa 



1 - 



(l0 - W 22 + Po7) + 



1 2 V 12 + 1^23^31 

3 ?7 



(86) 
(87) 



E\ _ Vyi _ % 1 JSi _ ^13 1 gg _ ^23 _. 1 

E 2 v 2 \ ' r\ E 3 v 3X A' £ 3 i/ 32 ' // 

The equations governing the equilibria of unshearable cylindrical shells are 
arrived at when the above constitutive equations are inserted into the balances 
(|6"5|) . In their general form, those equations are complicated to solve analytically; 
we choose not to list them here. However, certain highly symmetric problems 
admit simple and explicit solutions. We deal with such problems in the next 
section, while the noticeable simplifications obtained when the approximations 
judged appropriate for thin and slender shells will be introduced and exemplified 
in Section 5. 



4 Cylindrical Shells: Axisymmetric Boundary- 
Value Problems 

A boundary-value problem for a cylindrical shell is axisymmetric if the load and 
confinement data induce solution displacement fields whose physical components 
u<i> are all independent of the circumferential coordinate $, that is to say, in 
view of ASH), if 

m<i> = a<i> — (w' , u<2> = ( 1 + — | a<2>, u<3> = w + £7, (88) 

V Po) 

where a prime denotes differentiation with respect to x\, the only space variable 
from which all of the parameter fields a<i>,a<2>, and w, may depend. When 
the displacement field has the form (|88|) . 
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(i) the strain components (|63l) take the simpler form: 



E<u> = a<i> - C,w" , 
E<i2> = E<2i> = i ( 1 + — ] a<2>', 



E<22> = \ Po ( 1 
i?<33> = 7; 



Po 



Po 
-1 



(89) 



O + Ci), 



(ii) the constitutive equations 
become: 



for the force and moment components 



F<11> = 2 ±^1 



(1 - ^23^32) \p a<i>' - ^ £2w "^j 



Po A 

(v 2 l + f23^3l)w + (^31 + ^21^32)^07 



£ E 2 
F<22> = 2 



+ Po 



Po A 

V2\ + ^23^31 



(1 - ^132^31) ^ ( 1 - log 1 _ ^ 



1 1 + — 

Po7 + log i _ 1 w 

Po Po / 



V 



le 2 



F<2i> = 2eG ( 1 + -— ) a<2>', 

3 Po / 



F<33 > = 2±^ 
Po A 



a<i>' + (u 32 + ^31^12)^07 



f 31 + ^21^32 
A 



(1 - ^12^2l)p o 7 + 



p Q a<i> - yM) 



^32 + ^31^12 



■ W 



(90) 



and 

M<n> 



2e 3 Si 
3^ A" 



(1 - ^23^32)(Po w " - a<l>') - (1/31 + ^21^32)7 



4e 3 , 2e 3 . 

M<2i> = - — Ga<2> , M<i2> = - — Ga<2> , 

5 Po 3 Po 



M. 



<22> 



E 2 

2e i 



/ 1 i+f 

(1 - ^13^31) 1 - lo s 1 _ P J_ 

\ Po Po , 



1 . s 1 2^12 + ^23^31 „ 



(iii) the reaction- free equations ([55)1 reduce to the field equations: 
F<n>' + q a <i> = 0, 

(F<21> + p~ 1 M <2i>)' + q a <2> + J o~ 1 r Q <2> = 0, 
M<n>" — p Q x F<22> + q Q <3> + r Q <i>' = 0, 



(91) 



(92) 
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holding in the interval (— I, +1), plus the integral relation: 




[p a M <22> + F<33> — r a <3>) dxi — 



(93) 



(iv) the reaction-free boundary conditions consist in specifications at ±Z of one 
of the elements in each of the pairs (|7T|) 1 . ([7T]) 2 . and ([7T]) 1 : 



Insertion of (|M|) and (li?Tl) into (j9"2")l and yields the system of equations 
ruling axisymmetric boundary-value problems in our theory. We are going to 
solve this system for assignments of Neumann data corresponding, respectively, 
to problems of torsion, axial traction, pressure, and rim flexure. Interestingly, 
as we shall quickly demonstrate, this system splits into one equation for the 
circumferential displacement a<2> plus a system of three equations for the axial 
and radial displacements a<i> and w and the thickness stretch 7. 

Remark. As mentioned in closing Section ^. 51 once the equilibrium displacement 
field has been found, the axial distributions of reactive stress measures _F<3i> 
and Moi> consistent with boundary conditions compatible with (|7iP g and 
(|7TJ) 5 can be determined by a use of, respectively, the balance equation (|551) 1 
and the balance equation (|57l) , keeping into account (|93|) . The situations of our 
present interest occur when homogeneous Neumann data are prescribed at the 
boundary, so that the problems to solve are: 



M<3i>' = m(-l,+l), M<3i>(±0 = M<3i> = in (96) 

Note that equation (|95p i has the familiar structure of the moment balance in 
Bernoulli-Navier rod theory, with Foi> playing the role of the shear resultant, 
M<n> of the bending moment, and r a <i> of the diffused applied couples; both 
Foi> here and the shear resultant in that classic rod theory have reactive 
nature, as a consequence of one and the same unshearability constraint. 

4.1 Torsion 

From (|9"0)U and (|9"Tj) n, we have that 



(F<n>, a<i>), (i^<2i> + p D 1 M<2i>, a<2>), (A/<n>, u/); 



(94) 



F<3i> = M<n>' + r G <i> in (-£,+/), Foi>(±l) = 0, 



(95) 



and 




(97) 



with this, equation (I92IU takes the form of a second-order equation for a<2>: 




(98) 
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This equation accounts for whatever twisting about its axis an unshearable cylin- 
drical shell may have; it can be associated with boundary conditions specifying 
the values at x\ = ±1 of either a<2> or a<2>' . 

When the only applied load is a distribution of end tractions statically equiv- 
alent to two mutually balancing torques of magnitude 

T = (2np 2 )t, with t = O(s), (99) 

the thickness stretch 7 and the axial displacement a<i> vanish, and the 
only relevant equation, reduces to a<2>' = a constant, whose value is determined 
by the boundary condition: 

F<2i> + p~ x M <2i> = t = r T a<2>', r T := 2s ( 1 + ^\ G ; (100) 

it follows from (jlOOj) that a twisted shell of the type we study undergoes a 
rotation per unit length 

e:=p- 1 a<2>' (101) 
proportional to (/? rT) -1 . Moreover, given that the function a<2> must be odd, 

1 e~ 1 t 

a<2>{x 1 ) = - -5— — 11. (102) 

1 + % ZLr 

pi 



Remark. In line with (I75|) 2 . set 



q a <2> = q" l <2> = -S [ 1 + -— ] d<2> 



1^ 



in equation (|98|) . so that it takes the form of the classical wave equation: 



d 2 a<2>(x 1 ,t) 2 d 2 a<2>( Xl ,t) 2 C 1 

— w c — dx\ — = ' Wlth c= ^TT|f (103) 



Then, 

a<2>(xi,t) — ip(x\ + ct) + ip(xi — ct), 

and two twist waves propagate along the axis with speed |c|, the one in the 
positive direction the other in the negative direction. 

Remark. In all boundary-value problems we shall solve next - we recall, axial 
traction, uniform pressure, and rim flexure the twisting loads arc null. As 
a consequence, in all three cases, the second equation of system (|92|) and the 
boundary equations: 

(F<2i> + p- l M <2i>) (±i) = 

together imply that the construct (F<2i> + p~ 1 M<2\>) is identically null in 
[—1, +1]. Hence, by the constitutive relations ([90|) 3 and (|9T|)n. the circumferential 
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displacement a<2> has to have constant value, as is the case for a rotation about 
the x\— axis; we note that this, because of (l91~1) n. implies that 

M<2i>(±0 = 0, 

and we take a<2> = 0. In fact, as is customary in elasticity with Neumann 
data, in all three cases we expect to arrive at a displacement solution being 
unique to within an ignorable rigid motion. Moreover, given the common built- 
in symmetries, there will be no loss of generality in searching for solutions with 
a<i> an odd function of X\ € [—1, +1], and w even. 



4.2 Traction 

Let us take all distance forces and couples null and all boundary conditions of 
Neumann type and homogeneous, except for a distribution of tractions at the 
ends equivalent to two mutually balancing axial forces of magnitude 

P = (2np )p, with p = 0(e). (104) 

Two of the balance equations (|92|) are in force: 

F<h>' = 0, 

M<n>" - p~ 1 F<22> = 0; 

the accompanying boundary equations are: 

F<ii>(±0=p, M<n>(±0=0, M<ii>'(±2) = (106) 

the last one following from (p?5j) . Now, (11051^ and (I106p i imply that the mem- 
brane force F<n> is constant and equal to p in the closed interval [—1, +1]; then, 
making use of (1901^ . we obtain that the differential relation: 



(1 - ^23^32) (p a<i>' - ^e 2 w"^j + {v 2 \ + ^23^3i)w+ 



e-*p (1 ° 7) 

+ (^31 + ^21^32)Po7 = Po A 

must hold in [— Moreover, with the use of (|9T|) 1 and under the parity 
assumptions we made for a<i> and w, we see that the boundary conditions 
(I106IU take the common form: 

(1 - vaausa)(a<i>'(l) - p o w"(0) + (fai + ^21^32)7 = 0. (108) 

Finally, with the use of (|9T|) 1 again and of (J90j) 2 , and under the provisional 
assumption that the traction problem admits a solution with 7 a constant, 
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(|105[) 2 becomes: 



1 2 1 - ^23^32 ( n 2„„ll 



(fiw" - p o<i>') + 



3 ?7 

+ (1 - f 13*31) ^ f 1 " log Y^i) ^7 + ^ + ( 109 ) 

. V 2 \ + ^23^31 / . / | \ n 

H Po a <!> + (^32 + vz\vn)p l = . 

On eliminating a<i>' by means of (|10T[) . (|109[) yields: 
e 2 Po(l-^) »"" + £ 2 a/ + 6 w + c/)o A^+rfp o7 = 0, (110) 
where 

2(2/21 + 1^23^31) 
a := — , 

t — ^23^32 

; 3 / 1 1 + 1 ("21 +^3i) 2 \ 

:= i] 1 - z/i3^3i) ^T" log ^ i : , 

1 - ^23^32 V 2^- 1 - j- 1 - f 23*32 y 

_ 3(^21 + *23*3i) 

^ (l-^23^32) 2 ' C 111 ) 

3 / / 1 1 + 

d := 77(1 - 1/131/31) 1 - ttt" log " 



l-*23*32 ^ y 2 j- l ~jr o/ 

(*31 + *21*32)(*21 + *23*3l) , , 

+ ?/(//32 + *31*12j 



1 — ^23*32 / 

The general solution of the homogeneous equation associated with (|110[) has the 
form: 

Wh{x\) = ci(exp(aixi) + exp(— a\X\j) +c 2 (exp(a 2 a;i) + exp(— a 2 xi)) , (112) 
where 



a+Ja?- Ab§ (l - 1 g) -a - Ja* - 4 & f| (l - Jg) 



2 V \ ' / 9 

a a := ! t — , a 2 := 



(113) 

With this, we write: 

(c s~^~p d \ 
6 A "2^i" + 6 7 )' ^ 114 ' ) 

with w p the constant solution of (IllOp . 
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With a view to determining the coefficients c\ and C2, we firstly return to 
the boundary condition (I106|) 3 . that, when combined with (f9T|) . reads: 

p w"'(l) — a<i>"(l) = 0; 

in addition, by differentiating (|107|) and invoking continuity up to the boundary 
of the resultant expression, we obtain that 

ll fl \ 1 2 mfi\ ^21+^23^31 ,, n 

p a a<i> {I) = -e w {I) w (I); 

O i- — ^23^32 

the last two relations together imply that 

1 l£ 2 \ ,il fl\ i ^21 + ^23^31 1 i rn n / 11( -s 

1_ o-2 }Po w (0 + -j w(0=0. (115) 

3pg/ 1- ^23^32 Po 



Secondly, on eliminating a<\>'(l) in (I108|) by means of (|107p , we obtain: 

1 l£2 \ II (1\ , ^21 +^23^31 1 m 1 a n filfil 

3 Po/ 1 - ^23^32 Po (1 - ^23^32) 2E\ 

On taking (|114l) into account, the system of equations (|115l) and (|116|) deter- 
mines the coefficients c Q in (| 1 12p : 

exp(2ai I) — 1 

Cl = Ki 



C2 = «2 



ai ( exp(2ai exp(2a 2 i) + l) - a 2 ( exp(2ai + exp(2a 2 - l) ' 

exp(2a2 I) — 1 



a 2 ( exp(2a 2 - l) ( exp(2«i + - at ( exp(2a2 Z) + l) ( exp(2ai Z) - l) ' 

(117) 

with 



K-2 



\b \ 1-^23^32 / 2El l-V 



1 r° \ 3 pi J 1— 1/33V33 



ft K M ^21^23^31 ) A £ 1 P I "21+"23"21 d 

"1^6 V 1 1-1^23^32 ) 2Ex 1-1/23^32 6 ' 

™ 2 „ fix lE 2 "l i/l" 1 ^21+^23^31 

«2 p o ^l+3^J+Po l^^J 



(118) 



note that both the coefficients c Q depend on 7. 

Having found the form of the radial displacement w in [—1, +/], we revert to 
equation (I10T|) to find the axial displacement a<i>. A simply calculation yields: 



a<i>(xi) = 



! ^21 + ^23^31 £ \ ^ £ V 



1-^23^32 6/ 2i?i 



-( 1 - /'-:;/' ;.. ) ' ( t(^21 + ^23^3l) ~ (^31 + ^21^32) j 



1 



3p a aia 2 



12 2 ^21 + ^23^31 \ i i x / u 

«2Ci \e a.\ —6 ( exp(aixi) — exp( — a,\X\ ) )- 

' 1 - ^23^32 ' 



12 2 o "21 + ^ 23^31 \ I i , , ,\ 

«ic 2 e a 2 —6 — ( exp(a 2 a;i ) — exp(— a 2 xi j ) 

1 — ^23^32 ' 
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(needless to say, the parity condition a<i>(0) = is satisfied). The one task 
remaining is to find the constant 7. This we do by a sequence of manipulations 
of the integral balance equation (|93l) : 
(i) using and we give equation the form: 



(w + p j) - -e 



1 2 ^12 + ^23^31 // 



, n ^ , 1/31 + ^21^32 / , 1 2 /A , / \ I n 

+ Ml 1 - V\2V2\)P -1 ~\ I p a<i> - -£ W I + (^32 - 1^31^12 )w 1=0; 

(119) 

(ii) on expunging the construct (p D a<i>' — ^e 2 w") by means of (jTUTJ) , we have: 
2e 2 a w' h (l) + ai / w h =2l (p Q ( A^— - + a 2 7) - aiu^ , 



where 

1 1 - ^23^31 



a 



3 (V3I + ^21^32)(^12 + ^23^3l) 



1 - ^23^32 ( . H , ( 1 1 , 1 + p 

ai := 7/ ■ !/ 2 i + ^32^31 - ^23 + ^31^12 - (1 - f 13^31) 1 - 7TT 1o S 



^31+^21^32 V " \ 2f- l ~t~ ) 

a 2 ■— rjA — ^3^32 [ (1 _ j, 13i , 31 ) [ 1 - J_log- p -f- ] +/x(l - i/ 12 i/ 21 ) ) + 

^31 + ^21^32 y y 2^- 1 - J 

- (^31 + ^21^32) ; 

(iii) in view of (|114[) . we end up with: 

+ - — fai _1 ci(ai + e 2 a ai 2 ) sinh(aiZ) + a2 _1 e 2 (ai + e 2 a a 2 2 ) sinh(a 2 /) ) , 

lPo V 

(120) 

an implicit equation for 7. Figures [6][8] help visualizing some relevant features 
of the analytic solution we just constructed. 
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M<n> 




Figure 6: Axial Traction Problem: qualitative diagrams of bending moment 
M<n> and shear force F<3i>. 



f 


f ^ \w( Xl )-w(0)\ 













Figure 7: Axial Traction Problem: radial displacement w. 



Figure 8: Axial Traction Problem: cartoon visualization of deformed and unde- 
formed shapes. 
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4.3 Pressure 

We now let the cylindrical shell we study be subject to a uniform pressure 
g <3> = vj = 0(e), all the other applied loads being null. Accordingly, the field 
and boundary equations (|105l) and (I106P are replaced by respectively, 

F<n>' = 0, 

M <n>" - p~ l F <22> + zu = 0; 

and 

F<n>(±0 = 0, M<n>(±Z) = 0, M<u>'(±Z) = 0. (122) 

Equations (|12ip i and p22p i imply that F<n> is identically null, whence, in 
view of $U§ 1 , that 

/ 1 2 II v 21 + ^23^31 ^31 + I / 2lI / 32 r i , ;i 

p a<i> - -e w = — w p„7 over 

3 1 - ^23^32 1 - ^23^32 

(123) 

With this, d9TJ) x and ((90]) 2 , equation (|12ip 2 takes the form: 

e 2 p 2 (l-~^ w"" + e 2 aw" + bw-cp 2 o A^^ + dp o7 = 0, (124) 



where 



1 — ^23^32 

the other constants being defined by (|1 1 1[) x 2 r The general solution of the 
homogeneous equation associated to (|124|) has the form (|112|) : 

Wh{x\) = ci(exp(aixi) + exp(— aix{)) + c 2 (exp(a 2 x 1 ) + exp(— a 2 x{j) , (125) 

with the coefficients as given by (|113[k then, 

i \ ( \ ( c e~~ x vj d \ . . 

w(zi) = Wft(xi) + ui p , w p := p Q I -p c A - -7 I , (126) 



with w p the constant solution of (|124|) . The boundary conditions (|122|) 2 3 be- 

1 l£2 \ "fl\ , ^21 +^23^31 1 m n / 10 „x 

1 - o^r P„» (0 + ~, w(l)=0, (127) 



1 l£ ^ "'m i ^21 + ^23^31 1 / m n , ,. WM 

l--—\p w (Z) + - w(Z) = 0, (128) 

3 Po/ 1-^23^32 P 

and are expedient to determine the constants c\ and c 2 : 

exp(2ai I) — 1 



Ci = Ki 



C2 = «2 



a x ( exp(2ai I) - l) ( exp(2a 2 + l) - a 2 ( exp(2ai + l) ( exp(2a 2 I) - l) ' 

exp(2a 2 i) — 1 

a 2 ( exp(2a 2 Z) — l) ( exp(2ai I) + l) — ai ( exp(2a 2 Z) + l) ( exp(2«i Z) — l) ' 

(129) 
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with 

[CM 1^21^23^31 ) „ ASZ—E. — "21+"23"21 ) 

! \f> \ 1-1^231^32 y ^° 2£i l-y23f32 b ' ^ 



«2 I 

Ki := 



ai 2 p Q (1 + ±4 J + Po 1 ^ 1 ^ 23 ^ 31 

x ro V 3 / 1 — ^23^32 



«2 



rvi I £ M — i y 2H / 23i^3i ) „ A e 1 ro _ 1^21 + ^23^21 rf-, | 

^fe \ 1-1^23^32 y Fo 2Bi l-y23f32 



(130) 



a 2 2 P (l + ig)+Po 1 ^ Sf 



Finally, by a sequence of steps completely analogous to the one leading to (|120[) , 
the integral balance (|93|) yields an implicit equation for the constant 7: 



ai - + a 2 ) 7 =ai-p D A 
2 



6 ro 2Ei 

+ - — ^ai _1 ci(ai + e 2 a ai 2 ) sinh(o!i/) + a,2^ x ci{a\ + e 2 a a2 2 ) sinh(a2/)^ . 

(131) 



M<a> 



F<31> 





Figure 9: Pressure Problem: qualitative diagrams of bending moment M<u> 
and shear force Foi>. 



1 w(xi) 










-I 


1 



Figure 10: Pressure Problem: radial displacement w. 
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4.4 Rim flexure 

Lastly, we consider the case when a uniform distribution of bending couples 
m = 0(e) per unit length is applied at both rims of the cylinder. This time, 
the field equations we have to satisfy are: 

F<n>' = 0, 

( 132 ) 

M<u>" - p^F <22> = 0, 

with the boundary conditions: 

F<n>(±/) = 0, M<n>(±l)=m, M'<n>(±Z) = 0. (133) 

Just as in the pressure case, equations (|132j) x and (|133p i imply that (| 123[) holds; 
moreover, with the use of (jOTjl -p (p?!]]^ . and (II 23[) . equation (ll32P o takes the form: 



: 2 /r ( I - ™ ) -! : 2 o ,r" + hw i dp,,-. = ! 134 ) 



(cf. (1124|) ). where the constants are defined in (|lllj) . We set: 

W(xi) = - ^poJ + Whixi), 

with 

Wh(xi) = ci(exp(aixi) + exp(— a\X\)) + c 2 ( exp(a 2 xi) + exp(— a 2 xi)) , (135) 

the constants a\,a% being given by (|113p . On writing the boundary conditions 
(1133jU 3 in terms of displacements: 



l ^ \ ,i n \ , ^21+^23^31 1 m 3 A m Po 

3 Po / 1 - ^23^32 p 2 £ J i/i 



1 l£ 2 \ Wm . "21 + "23^31 1 / m n / 1Q „s 

1_ q^2 ^ (0 + -j w(0=0, 137) 

3 Pi J 1~^23^32 Po 

we determine the constants ci and c 2 



. _ f_i /3 mp Q d \ a 2 (exp(2a 2 Z) - 1)) exp(aiZ) 

Cl _ 1 \2 e 3 ^ ~ b 1 ) ai(exp(2a 2 + l)(exp(2aiZ) - l) -a 2 (exp(2a 1 + l)(exp(2a 2 Z) - l) ' 

„ _ t_i ( 3 mp d \ ai(exp(2o;iZ) - 1)) exp(a 2 Z) 

C2 ~ 2 V2 e^" _ 6 7 a 2 ( exp(2aiZ) + l) ( exp(2a 2 Z) - l) - ai ( exp(2a 2 Z) + l) ( exp(2aiZ) - l) ' 
where 

% (, 1 e 2 \ 2,-1 ^21 + "23^31 

fc l = U-q— I PO^ +P "I . 

V 3 pi J 1 - f 23^32 

f A 1 £ 2 A 2,-1 "21 + "23^31 

V 3 pi J 1 - f 23^32 
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Once again, the integral condition balance yields an implicit equation for 
d \ 2 



axT + a-2 |7 = 



] 7 = - — (at\ 1 ci (ai + e 2 a ai 2 ) sinh(aiZ) + «2 1 02(0.1 + e 2 a a2 2 ) sinh(a20 

/ l Po V 

(138) 



42 



M<n> 



F<31> 




II 
1) 



Figure 11: Rim Flexure Problem: qualitative diagrams of bending moment 
M<n> and shear force F<3i>. 
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Figure 12: Rim Flexure Problem: radial displacement w. 



Figure 13: Rim Flexure Problem: cartoon visualization of deformed and unde- 
formed shapes. 
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5 Cylindrical Shells: Thinness, Slenderness, Con- 
traction Moduli, and Stiffnesses 



As all figures from ([5]) to (I13[) make evident, a phenomenon of boundary lo- 
calization takes place, whatever the shell's thickness, in the boundary-value 
problems of traction, pressure and rim-flexure we have solved analytically. This 
phenomenon is more and more pronounced as the shell's length grows ceteris 
paribus. In our opinion, this fact amply justifies the use of the much simpler 
formulas for slender shells that we derive in this section. 



5.1 Thin shells and slender shells 

We term thin a cylindrical shell of diameter 2(p a + e) and length 21 if s/p a <C 1, 
slender if p a /l <C 1; most of times, p a is indeed smaller than /, so that a thin 
shell is slender as well (and e/l -C 1). For p a /l fixed, a thinner and thinner 
shell looses its bending and twisting stiffness more and more; in the limit for 
e/p a — > 0, it responds to loads as a tubular membrane. On the other hand, 
for both e and p Q fixed, longer and longer shells become slender and slender, 
without loosing their shell-like response. 

To exemplify the simplifications ensuing from taking large-thinness limits, 
we take the integral balance (|93|) . that has the same displacement form (| 1 1 9[) in 
all three boundary- value problems of traction, pressure, and rim flexure. Since 

/ lin Ln2^ l0g r 1 ^ = 1 ' (139) 
e/ Po ^0 2— 1 - — 



we see that (|I19j) reduces to 



p{l- v 12 i>2i)p "i + n x (^3i + ^21^32) -j-a<x>{l) + (i>32 - 2/31^12)^7 / w = 



21 



+1 



-/ 



where of course the values of both a<i>(l) and J_ ; w are problem-dependent. 

Remark. CNTs, no matter if single- or multi-wall, are as a rule slender. However, 
they may be thin or not, in the sense of the above definition [5{. Of course, all 
shell theories concern thin objects, but the thinness notions they are constructed 
upon may differ (see the discussions in [T2] and [13]); in particular, those notions 
need not be expressed in terms on one purely geometrical aspect ratio. Our shell 
theory works whatever that ratio, because its subtler and more complex notion 
of thinness is the one typical of the method of internal constraints, a method to 
derive the mathematical models of linear structure mechanics firstly sketched 
in0. 



5.2 Axisymmetric equilibria of slender shells 

Hereafter, we display the slenderness approximations of the solutions to the 
fundamental problems analyzed in the previous section, except for the torsion 
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problem, where no such approximation is in order, because the solution does 
not depend on the cylinder's length. 

(i) Traction . Looking at (|117p . it easy to conclude that 

lim Ci = lim ci = 0; 

l— >oo I— ► oo 

then, (|112j) takes the form: 

w h {xi) = 0. 

Moreover, (|120p yields that: 

7 = "7^, 7»^rA. (140) 
t&i a,2 + ai ^ 



With this, one hnds that 



e 1 p _ ( c 1 + an 

-T^T, w P-= Po\j + — 



w ( Xl )^-u, p —^, tg p := Po (-+ _ , A. (141) 



To arrive to a large-slenderness approximation for the function xi i-> 
a<i>(xi), the remaining unknown of our problem, we turn to (|107j) . Un- 
der the present circumstances, that equation yields the value of the axial 
strain in a slender shell: 



A := a<i>'(a;i) = 



_ A + p o 1 W p (v 2 l + ^23^3l) + 7(^31 + ^21^32) g 1 P 



1 — ^23^32 2£i 

(142) 



(ii) Pressure . This time, we look at (|129j) to conclude that, in limit for / — > 00, 
both constants c a tend to zero. Thus, once again, w(xi) = w p , and we 
find that 



7 = 7^, 7:=^A, (143) 



and 



^^Wj,^—, w p --=Po ( 144 ) 

Finally, a use of f|123[) yields: 

, ( s_ P c 7 1 Wp(^21 + ^23^3l) + 7(^31^21^32) PoS^rU 

a<i> (xi) = — ; — . (145) 

t — ^23^23 ^1 

(iii) Rim Flexure . Given that lim t\ = lim cq, = , it easy to show that the 
displacement field is everywhere null, in the large-slenderness limit. 



45 



5.3 Effective contraction moduli 

Let us now regard a slender cylindrical shell as a three-dimensional rod-like 
body, in short, a probe. On defining the cross-section strain measure: 

E( Xl ):=±J +e a(C)E(x u Q, (146) 

we have from ([55)) that 

E<n> — a<x> , E<i2> = E<2i> — i a<2>', E<22> = p~ l w, E<33> = 7. 

(147) 

Moreover, (|88|) 3 implies that the deformed external radius of a shell of unde- 
formed external radius r$ xt ~ p + e is r ext = p Q + w + £7, so that 

^ — ^ ~ T ° Xt = E<22> + 0(e). 



r ext 



Motivated by this observation, by an effective contraction modulus we mean: 

• in case of axial traction, 

VcA ._ g< 22 > _ Po 1 Wp{l - V2ZVZ2) ^ 4g ^ 

E<u> A + po 1 w p (v 2 i + ^23^31) + 7(^31 + ^21^32) ' 

where w p and 7 are given by, respectively, (|141[) and (|140l) (note that this 
definition is directly reminiscent of the laboratory procedure followed to 
measure Poisson's modulus for isotropic materials); 

• in case of uniform inner pressure, 

E<n> Po 1 w p {v 2 i + ^23^31) +7(^31 + ^21^32) ,,, n , 
v C p-=-= = -1- ,., : ; ( 14 9) 

E<22> p W p (l- V 2 ZV Z2 ) 

where w p and 7 are given by, respectively, (|144l) and (I143[) . 
Both vca and vcp depend on thickness, through, respectively, w p ,j and w p ,j. 



5.4 Effective stiffnesses 

It remains for us to introduce suitable notions of effective traction and torsion 
stiffnesses for a slender cylindrical shell regarded as a probe. This we do by 
mimicking the relative familiar formulas from one-dimensional rod theory. 

As to effective traction stiffness, given that the traction stiffness of a rod is 
defined to be (axial load) /(axial strain), we set: 

SA ■= t> ( 15 °) 
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whence 



(1 - 



ExA{e) , A{e) := 4np a e, 
(151) 



SA = 



A + p 1 w p {v 2 \ + v-izvzi) + 7(^31 + ^21^32) 



where we have made use of (|104[) . (|142l) . and (|147p i . and where A(e) is the area 
of the shell's cross-section (needless to say, sa depends on e also through w p 
and 7). 

As to effective torsion stiffnesses, we recall that, for a twisted rod, one takes 
it to be (torsion moment) /(twist per unit length); accordingly, we set: 



where J(s) approximates to within 0(e 2 ) terms the polar inertia moment of the 
cross section, and where x{ £ ) = 1 + 0(e 2 ) is a sort of torsion factor F°l 

6 Cylindrical Kirchhoff— Love Shells 

Recall from Section 12.21 the Kirchhoff-Love representation (fT4"]) for the displace- 
ment field: 

u kl(x,C) = s A(xX) a (x) + w(x)n(x) — ( s Vw(x), a(x) ■ n(x) = 0, 

and compare it with the more general representation (111[) : in both cases the 
unshearability constraint is imposed implicitly, together with, but only in the 
first case, inextensibility of the fibers orthogonal to the middle surface. Thus, 
the thickness of these shells does not change whatever the applied loads, making 
the simpler Kirchhoff-Love theory suitable, in our opinion, for application to 
single-wall CNTs. 

Formally, to recover the representation (fl4"|) . one only has to take 7 = in 
(llip . To adjourn the developments of Sections 3, 4, and 5, this measure has to 
be accompanied by a few other adjustments that we now categorize and detail. 

(i) (Balance Assumptions) The variations entering the Principle of Virtual 
Powers must have the following form, to be compared with (|64|) : 



Consequently, the integral balance equations (1551) 1 disappears, and the 
same happens with the last power-conjugate pair in each of the boundary 
conditions ([71]) and ([M]) . 

Note that l|153|l holds whatever the slenderness of the shell under consideration. 




(152) 



(153) 



s v = V1C1 + v 2 n' + v 3 n, 

V = -1> 3 ,1 Ci + Po 1 ^ - v 3 , 2 )n', 



(154) 
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(ii) (Constitutive Assumptions) The constraint space (j4"Tj) reduces to 

9 = span( W h i = 1,2,3); 

<S'<33> is then reactive and only four constitutive moduli survive. It is not 
difficult to show that ([50)1 and (|5"Tj) must be replaced, respectively, by 



= A ( Ei Wi <g) Wi + E 2 W 2 ® W 2 + 2A G W 3 ® W 3 

(155) 

Eii/ 2 i(Wi ®W 2 +W 2 ® Wi) 



and 

A := 1 - 1/12 ^2i- (156) 

This result is achieved by the first of various applications to follow of a 
procedure that consists in taking the limits for £"3, A, and fx, tending to 
infinity, and in setting v a3 — v 3a = (a = 1,2). It is important to 
realize that, due to the first of ([57)) that we here repeat for the reader's 
convenience: 

T = (157) 

Hi2 V 2 \ 

(I155|) and (|156|) integrate a 4-parameter constitutive representation. 



(hi) (Torsion Problem) The solution given in Section 14.11 does not change, be- 
cause it does not involve any of the constitutive moduli that take singular 
values in the case of Kirchhoff-Love shells. 

(iv) (Traction Problem) The two governing equations are, mutatis mutandis, 
(I107P and (|110[) : specifically, the equation corresponding to (|10T|) is: 

Po a<i>' - V + v 2 i w = Po A E -^- , (158) 

while the one corresponding to (II 10[) reads: 

e 2 p 2 (l - ~) w"" + e 2 aw" + b w + c Po A ^ = 0, (159) 

where the constants have the following expressions, that can be recovered 
from (fTTTj) : 

a:=2v 21 , b := 3rj I -^j- log _ ^ - ^12^21 ) , c:=iv 21 . (160) 
As to equation (|114p . it is replaced by 

w(x!) = w h (xi) + w p a, w p a = -p Q v\2 6(e) — — , (161) 
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where 

5(e) := — ^_ 2 ^ 21 ^ (note for later use that lim 5(e) = 1). 

5V log rf-^21 £/p °^° 

Po Po 

(162) 

For slender shells, w(xi) = w p and (|158l) allows to conclude that 

0<l>(Xi) = (l - ^12^21 (1 - 5(e)) ) "TT^ Xl. (163) 

(v) (Pressure Problem) The two relevant equations, (I123[) and (|124[) , become, 
respectively: 

P a<!>' - t^V' = -^21^, (164) 

and 

^(l-^) w""+e 2 aw" + bw-3p 2 o A £ -^=0, (165) 

where the constants a and b are the same as in equation (|160l) . Equation 
(|126[) is replaced by 

w(xi) = w h (xi) + w p p, w p p = p 2 S(e) £ CT . (166) 

By using (|164[) . we conclude that, for slender shells, 

o<i> = -i/ ai *( e )M^ ari . (167) 

(vi) (Effective Contraction Moduli. Effective Stiffnesses) As to contraction 
moduli, formulas (|148l) and (I149p become: 

5(e) 

vca = -r. 777-7 V12 and v CP = v 2 i- (168) 

1 - ^12^21 (1 - 5(e)) 

As to traction stiffness, (|151l) reduces to 

sa = z K 77^7 E 1 A(e); (169) 

1 - ^12^21 (1 - 5(e)) 

torsion stiffness continues to be given by formula (|153[) . 

(vii) (Determination of Effective Constitutive Moduli) In principle, a set of 
simple torsion, traction, and pressure, experiments or simulations allows 
to deduce the values of st,sa, vca, and ^cp, from the prescribed values 
of applied torque, axial load, and inner pressure and the measured or 
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computed values of 0,A, and T. In the large-thinness limit, one would 
have: 

GJ(e) = s T , E 1 A(e)=s A , "12 = "ca, v- lx = v C p- l 

so that the values of the four constitutive moduli G, E\, V12, and V21 would 
follow, provided one could measure or evaluate the values of the geomet- 
rical parameters p a and e. Now, given the well-known difficulties, com- 
monly referred to as the Yakobson's Paradox [16) . in choosing a repre- 
sentative value for the wall thickness of a SWCNT, we think it best to 
characterize the mechanical response of a Kirchhoff-Love shell by the con- 
traction moduli ^12,^21 and the effective constitutive moduli G := Ge, 
E a :=E a e, (a = 1,2). 

Remark. We find it appropriate to expand a little on this last point. For 
slender and thin cylindrical shells, we find that: 

- when subject to end traction, 

wt = -p v\2^^, a<i> T = (170) 
2E X 2E 1 

by way of equations (|161[) and (I163[) ; 

- when subject to inner pressure, 

2 m I Po W 

w P = p -^, a<i> P = -v 2 \—^, (171) 
2E2 2E2 

by way of equations (|166[) and (|167[) . 

Equations (|170[) and (|171l) can be regarded as a system of four equations 
in the unknowns E±, E2, V12, ^21, whose solution is: 

Z P p 2 w w t a<i>' p 

2a<i> T zfflp p Q a<i> T wp 

(172) 

Interestingly, the consistency condition (|87p i . which can now be written 

as _ 

-Bi _ ^12 
E 2 ^21 ' 

implies that 

p to 



wt a<i>p 

a relation that can be used to check whether the present simplified form 
of our shell theory is applicable. 
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7 Conclusions and further developments 



We have given a detailed presentation of a theory of linearly elastic orthotropic 
shells with potential application to the continuous modeling of carbon nan- 
otubes. The novelty of this theory resides in two features: (1) the type of 
orthotropic response we have selected seems suitable, with minimal tuning, to 
capture chirality, not only in the extreme cases of zig-zag and armchair SWC- 
NTs, but also when it varies in an essentially undetectable manner from wall to 
wall of a MWCNT; (2) the possibility of accounting for overall thickness changes, 
that should be almost exclusively due to changes in inter-wall separation. As 
a matter of fact, the referential thickness of an ideal MWCNT is dictated by 
inter-wall forces of van der Waals type, whose action may be modeled essen- 
tially in two ways: either they can be regarded small with respect to applied 
loads, and ignored altogether; or they can be thought of as inducing a referential 
equilibrated stress state that should be taken into account when studying the 
effects of boundary conditions, no matter if hard or soft: our sounding of this 
latter approach is encouraging. 

In addition, we have proposed a simpler version of the theory, in which 
orthotropy is preserved but thickness changes are excluded in all admissible de- 
formational vicissitudes; we believe this simpler theory to fit SWCNTs, whose 
effective thickness when regarded as cylindrical shells should not change appre- 
ciably when loaded no matter how evaluated. 

Presuming a rather complex material response requires specification of a 
number of constitutive parameters - seven when thickness changes are allowed, 
four when they are not. Luckily, another feature of our present theory is that, 
in both its versions, it leads to a number of significant boundary- value problems 
that can be solved explicitly in closed form. These problems are: torsion, axial 
traction, uniform inner pressure, and rim flexure; were their solutions coupled 
with the corresponding measurements and/or simulation results, applicability 
of our theory could be unequivocally assessed and all constitutive parameters in 
it uniquely determined. It is not difficult to paste explicit solutions of the type 
we here derived for two or more coaxial shells, both in statical and dynamical 
situations; we are currently developing this line of research with a view to a 
better understanding of van der Waals forces and the role and evolution of 
defects [2]. 
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